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® Image processing method and device for performing that method. 



© An image processing method is provided ac- 
cording to the invention to remove areas from an 
image made by x-irradiation, said areas correspond- 
ing to absorption fitters in the x-ray beam. Because 
such filters, also known as x-ray shutters, are in- 
serted into the x-ray beam from the; outside of the x- 
ray beam, x-ray shutter-areas extend from the pe- 
riphery of the image. Therefore, once the edge of an 
x-ray shutter-areas has beep determined, the part of 
the image located between the periphery and the x- 
ray shutter-edge is removed. X-ray shutter-edges 
are determined by locating transition-points in the 
image having brightness-transitions having maxima 



of both the first and second spatial derivatives of the 
brightness-value. An x-ray shutter-edge is then de- 
termined by fitting a curve through the transition- 
points. The curve fitting Is performed with the use of 
weight-factors that a re representative for the likeli- 
hood of an transition-point to be on an actual x-ray 
shutter-edge. This likelihood is larger for an transi- 
tion-points which is closer to the periphery of the 
image, for transition-points having low image rough- , 
ness towards the periphery of the image. Further- 
more, an transition-point having a large height of its 
brightness-transition is more likely to be on an actual 
x-ray shutter-edge. 
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The invention relates to an image processing 
method formed by x-ray imaging, and containing a 
portion corresponding to an x-ray absorption filter, 
the invention also relates to an image processor to 
process an image formed by x-ray imaging, which 
contains a portion corresponding to an x-ray ab- 
sorption filter. The invention further relates to an x- 
ray examination apparatus comprising an x-ray ab- 
sorption filter between an x-ray detector and an x- 
ray source for irradiating an object to form an x-ray 
image on the x-ray detector. 

An image processing method of said kind is 
known from the US patent US 4 670 896; therein 
are disclosed an x-ray examination apparatus and 
the use thereof to form on the x-ray detector the x- 
ray image containing a portion corresponding to 
the x-ray absorption filter. 

The cited reference teaches that a substantial 
difference in image-contrast can occur between 
parts of an x-ray image formed due to irradiation of 
the object by x-rays. Brightness differences for- 
ming image contrast arise e.g. in peripheral angiog- 
raphy because part of the x-ray beam does not 
penetrate the object, notably a patient to be exam- 
ined, and is incident directly on the x-ray detector, 
or for example in cardio-examination ^ because the 
object exhibits substantial differences in x-ray ab- 
sorption, viz. a patient's heart having comparatively 
high x-ray absorption is surrounded by lung-tissue 
which is relatively transparent for x-rays. In order to 
make the dynamic range of the x-ray image co- 
incide as well as possible with the contrast in the 
region of interest, x-ray absorption filters e.g. hav- 
ing the form of x-ray shutters, are arranged be- 
tween the x-ray source and the object The x-ray 
shutters are positioned to intercept a part of the x- 
ray beam so as to shield relatively transparent 
parts of the patients anatomy form the x-rays. The 
use of such x-ray shutters leads to diminishing an 
x-ray doses to which a patient is exposed to pro- 
duce an x-ray image. Moreover, an x-ray dose 
whereto a person who operates the x-ray examina- 
tion apparatus is unintentionally exposed is dimin- 
ished by the use of x-ray shutters. In the x-ray 
image x-ray shutters appear as regions having ex- 
tremely high absorption. Consequently, in inverted 
images, as are often used by radiologists for ex- 
amination, x-ray shutters appear as bright areas 
which distract attention from medical information. 
Furthermore, when a plurality of images made by 
x-ray irradiation is brought, together onto a hard 
copy, such as an x-ray film, a substantial portion of 
the area of the x-ray film is wasted when said 
images contain areas corresponding to x-ray shut- 
ters. Thus, making images, in accordance with the 
cited reference, by irradiating an object by an x-ray 
beam and employing the x-ray absorption filter 
arranged between the x-ray source and the x-ray 



detector has several shortcomings and disadvan- 
tages. ■ ■ / 

Therefore, it is an object of the invention to 
provide an image processing method for reducing 

5 harmful effects of areas in the image corresponding 
to the x-ray absorption filter on the quality of the 
image made by x-irradiation. 

This object is achieved by the image process- 
ing method according to the invention to process 

w an image formed by x-ray imaging, the image 
containing a portion corresponding to an x-ray ab- 
sorption filter, characterized in that the image pro- 
cessing method comprises modifying the image by 
selectively processing said portion. 

75 A processed image that does not show harmful 

effects such as distracting features in the portion 
corresponding to the x-ray absorption filter appear- 
ing in the image as x-ray shutter-areas is obtained 
in accordance with the image processing method 

20 of the invention since portions of the x-ray image 
corresponding to an x-ray absorbtion filter i.e. for 
instance x-ray shutter-areas are modified for exam- 
ple in that these are removed from the image. 
Then in the processed image the portion corre- 

25 sponding to an x-ray shutter area is absent. Pixel- 
values in the x-ray shutter-area may also be set to 
a predetermined grey-level or are given a non- 
distracting uniform colour-value. Therefore, in the 
processed image the portion that corresponds to x- 

30 ray shutter-areas is presented as an area having a 
non-distracting grey-tone or inconspicuous colour 
and consequently the processed image only con- 
tains image information that is relevant. 

It is a further object of the invention to provide 

35 an image processing method for identifying the 
portion corresponding to an x-ray absorption filter 
before selectively processing the image portion 
which corresponds to the x-ray absorption filter. 
A preferred implementation of an image pro- 

40 cessing method in accordance with the invention is 
characterised in that the processing comprises de- 
termining an edge of the portion for selecting said 
portion as a part of the image between the edge 
and a periphery of the image. 

45 Sections of the x-ray beam that are absorbed 

by the x-ray absorption filter may have a com- 
plicated shape corresponding to complicated ana- 
tomical structures of a patient to be examined. 
Therefore it is difficult and often quite impossible to 

so know beforehand or predict the location of an area 
which forms the portion of the image correspond- 
ing to said filter, i.e. an area of the irnage in which 
x-rays were absorbed due to the x-ray shutters. In 
accordance with the invention an image processing 

55 method is provided that can be performed indepen- 
dently of the x-ray examination apparatus and/or 
adjustments thereof that is employed to irradiate 
the object. Furthermore, it is often desired to pro- 
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cess the image separately from the x-ray examina- 
tion. Any information regarding x-ray shutter-posi- 
tions and/or adjustment will then be lost 

In images made by x-irradiation of an object 
one or several x-ray shutters are employed for 
avoiding overexposure in parts having low absorp- 
tion and being irrelevant from a medical point of 
view. Such x-ray shutters are without exception 
inserted into the x-ray beam from the edge of the 
x-ray beam towards the centre thereof. Conse- 
quently, when the (ocation of the edge of an x-ray 
shutter is determined it is certain that the entire 
region of the image from said edge towards the 
periphery of the image is an x-ray shutter-area, I.e. 
an area in the image where x-radiation has been 
absorbed by the absorption filters. 

When locations of edges of x-riay shutter-areas 
have been determined, a processed image wherein 
x-ray shutter-areas are reduced in that said x-ray 
shutter-areas are converted into a simpler form by 
deleting image information of x-ray shutter-areas or 
by setting pixel-values in areas peripheral to the xt 
ray shutter-edges to a predetermined grey-level or 
a non-distracting colour-value. A processed image 
that does not have distracting features in x-ray 
shutter-areas is obtained in accordance with the 
image processing method of the invention since x- 
ray shutter-areas are removed from the image 
and/or pixel-values in the shutter-area are set to a 
predetermined grey-level or colour-value. There- 
fore, in the processed image regions that cor- 
respond to x-ray shutter-areas are presented as 
parts having a non-distracting grey-tone or incon- 
spicuous colour and consequently the processed 
image only contains image information that is rel- 
evant. 

A preferred implementation of an image pro- 
cessing method in accordance with the invention is 
characterised in that the processing comprises re- 
moving at least a part of said portion from the 
image. 

A particular advantage of removing portions of 
a portion of forming an x-ray shutter-area is 
achieved when a collection of processed images 
are brought together on a hard-copy, e.g. an x-ray 
film. As irrelevant x-ray shutter-areas are discard- 
ed, less space on the x-ray film is required, so that 
on the same area of the x-ray film more images 
can be collected. 

A further preferred implementation of an image 
processing method according to the invention is 
characterized in that the processing comprises as- 
signing substantially uniform values to pixels in the 
portion in the image. 

Distracting portions in the image are avoided in 
that x-ray shutter-areas which do not contain medi- 
cal relevant information are converted into portions 
having a uniform grey-tone or an inconspicuous 



colour. Then the processed image is suitable to be 
considered by a physician who may comfortably 
focus his/her attention to the medical information 
contained in the image without being distracted by 

s the x-ray shutter-areas. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the detection of said 
edge comprises 

io a grouping step to select a group of transition- 
points in the image having a higher brightness 
towards the centre of the image and a lower bright- 
ness towards the periphery of the image, and a 
curve-fitting step to fit a curve through said group 

;5 of transition-points, said curve forming a repre- 
sentation of said edge. 

Transition-points in the image that lie on an 
edge of an x-ray shutter-area have high brightness 
towards the centre of the image and low brightness 

20 towards the periphery of the image. For each edge 
in the image a number of transition-points are 
determined and grouped into groups of transition- 
points for respective x-ray shutter-edges. Repre- 
sentations of respective x-ray shutter-edges are 

25 obtained by performing a curve-fitting method for 
each of the groups. The result is one or a number 
of curves being representative for x-ray shutter- 
edges in the image. 

A further preferred implementation of an image 

30 processing method in accordance with the inven- 
tion is characterised in that the grouping step com- 
prises a region-selection step to select a region of 
interest extending from a periphery of the image 
and a division-step to divide said region of interest 

35 into sections and a weight-assignment step to as- 
sign weights to said transition-points, a weight for a 
transition-point depending on image features ar- 
ound said transition-point and a weight for a transi- 
tion-point depending decreasingly oh a distance 

40 between said transition-point and said periphery. 

In an image made by x-irradiation it is known 
beforehand that x-ray shutter-areas extend from the 
periphery of the image because the x-ray shutters 
are inserted in the x-ray beam during x-ray expo- 

45 sure from the outside of the x-ray beam. Con- 
sequently, it is effective to carry-out grouping of 
transition-points in accordance with known topology 
of the image. A group of transition-points being 
representative of an x-ray shutter-edge h<as transi- 

50 tion-points located in a region of interest that ex- 
tends from the periphery of the image. A bright- 
ness-transition is a necessary but not a sufficient 
feature of a point in the image for being on an x- 
ray shutter-edge. Detection of transition-points is 

55 carried-out by dividing a region-bf-interest into sec- 
tions, notably elongated strips. Subsequently, 
brightness-values are averaged in a direction per- 
pendicular to the longitudinal axis of elongated 
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strips. Then, one or several points having a bright- 
ness transition are determined in each of the sec- 
tions. In each of the sections^ the brightness-transi- 
tion being closest to the periphery of the image is 
the most likely candidate for being on an x-ray 5 
shutter-edge/ This is implemented by assigning 
weight-factors to brightness-transitions found in 
each of the sections, each weight-factors depend- 
ing decreasingly on the distance between a bright- 
ness-transitions at issue and said periphery of the 10 
image. ; 

A further preferred implementation of an image 
processing method in accordance with the inven- ' 
tion is characterised in that the grouping step com- 
prises a first-estimate step to determine a first- is 
eistimate for a position of a transition-point in a 
section and an accurating step to determine a 
transition-point from said first-estimate. 

Processing speed is increased by performing a 
first estimate for a transition-point on the basis of a 20 
component of the image having only , low spatial 
frequency variations of brightness-values. Such a 
component is obtained by performing a local spa- 
tial averaging of brightness-values in the image. 
After a first estimate for a transition-point has been 25 
obtained. » the position of said transition-point can 
be determined more accurately by considering 
brightness-values of the image in a neighbourhood 
of the first estimate. Thus, it is avoided to a large 
extent to search in vain for transition-points in re- 30 
gions of the image where there are no substantial 
brightness variations. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the weight-assignment 35 
step comprises computation of image roughness, 
and computation of transition-height, and computa- 
tion of said weights depending on image roughness 
and transition height for respective transition-points. 

Brightness-transitions that relate to x-ray shut- 40 
ter-edges are distinguished from brightness-transi- 
tions that for example relate to anatomical struc- 
tures by taking image features around these trahsi- ' 
tipns into account Image roughness is defined as 
e.g. the statistical variance of brightness-values or 45 
e.g. the maximum of the differences of pixel-values 
of pixels having a predetermined distance in the 
image between them. It is found empirically that 
image-roughness is considerably lower in the low- 
brightness region in an image caused by an x-ray 50 
shutter. Also the height of a transition due to an x- 
ray shutter-edge is higher than a brightness-transi- 
tion due to anatomical structures. Therefore, selec- 
tion of brightness^transitions. relating to x-ray shut- 
ter-edges is made more effective if the weight- 55 
factors of the brightness-transitions are made de- 
pendent on the above image features. This depen- 
dence is such that a brightness-transition having 



image features indicating an x-ray shutter-edge are 
given higher values for its associated weight-factor. 
There can be several transition-points in a section, 
then transition-points located nearer to the periph- 
ery of the image add to image-roughness for tran- 
sition-points located nearer to the centre of the 
irnage. Consequently, a dependence which is de- 
creasing with the distance between the transition- 
points and the periphery of the image for weight- 
factors for transition-points within a section is ob- 
tained. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the weight-assignment 
step comprises computation of first and second 
spatial derivatives of brightness-values near re- 
spective transition-points and employing substantial 
concurrence of a maximum of said first derivative 
with a brightness-variation and with a maximum of 
said second derivative at the peripheral side of said 
brightness-variation for selecting the position of 
said brightness-variation as a transition-point and 
assigning said weight on the basis of said concur- 
rence. 

Brightness-transitions due to x-ray shutter 
edges have a particular image feature that distin- 
guishes x-ray shutter-edges from brightness transi- 
tions for example due to anatomical structures or 
vignetting. It is found that a brightness-transition of 
an x-ray shutter-edges has a maximum spatial first 
derivative of the brightness values near the bright- 
ness-transition in the image and has also a maxi- 
mum of the spatial second derivative of the bright- 
ness values slightly peripheral to the brightness- 
transition in the image. This feature is employed in 
selection of transition-points by assigning a low 
weight to transition-points not having a substantial 
value for their maxima of first derivatives and not 
having a maximum second spatial derivative slight- 
ly peripheral to the transition-point. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the weight-assignment 
step comprise? Computation of first spatial deriva- 
tives of brightness-values peripheral to respective 
transition-points, computation of respective average 
brightness-value of pixels towards the centre to 
said transition-points, computation of respective ra- 
tios of said respective first spatial derivatives and 
said respective, average brightness-values and as- 
signing said weights on the basis of said ratios. 

An x-ray shutter-edge gives rise to a bright- 
ness-transition in trie image stronger than contrast 
due to an anatomical structure. Another way of 
distinguishing an x-ray shutter-edge from anatomi- 
cal structures in the image consists of considering 
at a transition-point a gradient of brightness-values 
in a direction towards the centre of the image. A 
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value of such a gradient is subsequently connected 
to the average brightness-values between said tran- 
sition-point and a periphery of the image. If said 
gradient is comparatively strong and said average 
corresponds to high brightness then said transition- 
point is comparatively more likely to be on an x-ray 
shutter-edge. A weight for a transition-point is 
made dependent on the ratio of said gradient and 
said average, so that a transition-point is assigned 
a higher weight when it is likely to be on an x-ray 
shutter edge. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the weight assignment 
step is followed by a point-selection step for select- 
ing a transition-point having the highest weight in a 
section and retaining only said selected transition- 
points in said selected group. . 

Processing speed of the image processing 
method is increased by retaining as few transition-, 
points as possible. The transition-points that are 
retained should preferably be as likely as possible 
to be on an x-ray shutter-edge instead of being 
representative for anatomical structure. Such transi- 
tion- points are selected by sorting transition-points 
of a section on the basis of their respective weights 
and retaining in each section only the transition- 
point having the highest weight in its section. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised, in that said curve-fitting step is 
followed by a reassigning step for reassigning a 
value to a weight of a transition-point in depen- 
dence on a distance between said transition-point 
and a curve fitted through the group whereto said 
transition-point belongs* and in that the curve-fitting 
step is reiterated iemploying said re-assigned 
weight 

When anatomical structures are visible due to 
scatter radiation in an x-ray shutter-area, transition- 
points within the x-ray shutter-area may be se- 
lected in the point-selection step. The ensuing 
curve fitting step will then produce a curve to which 
selected points may have a comparatively large 
distance. This occurs because anatomical struc- 
tures have a more complicated structure than ap- 
proximately straight x-ray shutter-edges. To pre- 
vent that image structures that are e.g. due to 
scatter radiation are mistaken for x-ray shutter- 
edges, weight-factors of the transition-points are 
decreased on the basis of said distance. Transition- 
points which are close to the fitted curve have a 
higher likelihood to be actually on the edge of an x- 
ray shutter-area; therefore, weight-factors of transi- 
tion-points are adapted before performing a next 
iteration-step in correspondence with said likeli- 
hoods. The iteration is performed particularly effec- 
tively because an adequate data-set for starting the 



iteration is chosen as follows. For each section a 
peak-point is selected among the transition-points 
in said section. A peak-point is a transition-point 
that has the smallest average distance between 

s said peak-point and a predetermined number of 
neighbouring transition-points; empirically it has ap- 
peared that the choice of five neighbouring transi- 
tion-points is very adequate to determine a peak- 
point. The first estimate for a curve representing an 

10 x-ray shutter edge is then taken to be a curve that 
connects peak-points of respective sections. 

A further preferred implementation of an image 
processing method in accordance with the inven- 
tion is characterised in that the edge-determination 

is step is followed by an edge-verification step com- 
prising computation of a reliability-parameter for an 
edge being determined and comparison of said 
reliability-parameter with a predetermined threshold 
level. 

20 Even in the cases where no x-ray shutter was 
employed during irradiation the image processing 
method in accordance with the invention may pro- 
duce a curve allegedly identifying an x-ray shutter- 
edge. In such instances an anatomical structure 

25 would mistakenly regarded as an x-ray shutter. To 
avoid such errors, an edge-verification step is per- 
formed. Anatomic structures differ from x-ray shut- 
ter-edges in that the latter give rise to compara- 
tively straight curves. A plurality of said straight 

30 curves constitutes a, possibly complex, circum- 
ference of a section of the x-ray beam that is 
passed through the x-ray absorption filter. More- 
over, if a curve detected in a region of interest 
extending from one periphery of the image indeed 

35 represent an x-ray shutter edge, then most often 
another x-ray shutter-edge is present in a region of 
interest extending from an opposite periphery of 
the image. An average distance of a predetermined 
number of transition-points generates a quantitative 

40 measure of straightness of a curve. A reciprocal of 
said average distance is employed as a reliability 
parameter. If said reliability parameter has a value 
below a predetermined threshold-value then it - is 
decided that the detected curve does not for cer- 

45 tain relate for certain to an x-ray shutter-edge and 
the step of reducing the area peripheral to said 
curve consists of leaving the pixel-values in said 
area unaltered. If the reliability exceeds said 
threshold-value then a presence of an x-ray shut- 
so ter-edge in an opposite region-of-interest is veri- 
fied. If no such opposite x-ray shutter is found then 
the fitted curve is discarded and again the step of 
reducing the area peripheral to said curve consists 
of leaving the pixel-values in said area unaltered. 

55 A further preferred implementation of an image 

processing method according to the invention Is 
characterized in that for a plurality of transition- 
points said reliability parameter is computed and in 
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a situation where comparison of said reliability pa- 
rameter to said predetermined reference level in- 
dicates the presence of said portion corresponding 
to an x-ray absorption filter for less transition-points 
than a predetermined reference number the edge- 
verification step is supplemented by comparison of 
pixel-values of collections of pixels on different 
sides of the curve. 

When the number of transitions point retained 
for fitting the curve is less than a predetermined 
reference number then an ambiguous situation with 
regard to deciding on the presence of an x-ray 
shutter arises. Then according to the invention the 
average intensity of collection of pixels on either 
sides of the curve are compared. Should these 
average intensities differ more than a predeter- 
mined reference level than it is decided that an x- 
ray shutter is present because a high mean inten- 
sity on the side of the curve towards the centre of 
the image and a low mean intensity on the side of 
the curve toward the periphery of the image is 
compatible with the presence of an x-ray shutter. 
The predetermined reference number may be set 
to a very high value, for instance equal to the total 
number of pixels of the image, in order to achieve 
that supplementing the edge-verification step is 
performed in almost any situation. 

To reduce the computational burden said aver- 
age intensities may be calculated for comparatively 
small regions in the image, for example two vertical 
line segments an inner one with respect to the 
calculated curve lying towards the centre of the 
image and the other, outer one, lying towards the 
periphery of the image. 

To improve the reliability of the edge verifica- 
tion the maximum value of the second spatial 
derivative (Laplacian) of the pixehvalues on the 
outer segment is determined. If this maximum val- 
ue is found to be large, it indicates the presence of 
a anatomical structure towards the periphery of the 
accurate transition-point at issue and hence the 
likelihood of the actual presence of an x-ray shutter 
is reduced. 

A further preferred implementation of an image 
processing method according to the invention char- 
acterized in that plurality of transition-points said 
reliability parameter is computed and in a situation 
where comparison of said reliability parameter to 
said predetermined reference level indicates the 
presence of said portion corresponding to an x-ray 
absorption filter for less transition-points than a 
predetermined reference number the edge-verifica- 
tion step is supplemented by determining a mag- 
nitude of a first spatial derivatives of brightness- 
values of pixels in the image between the curve 
and the periphery of the image, said derivative 
being negative in the direction from the curve to 
the periphery of the image and comparing said 



magnitude to a predetermined reference value. 

When the number of transitions point retained 
for fitting the curve is less than a predetermined 
reference number then an ambiguous situation with 
6 regard to deciding on the presence of an x-ray 
shutter arises. Then according to the invention the 
gradient of the intensity for collection of pixels on 
either sides of the curve are compared to each 
other. As intensities, and thus pixel values, hardy 

10 ever decrease strongly within an x-ray shutter-area 
in a direction towards a periphery of the image, 
considering the most negative slope constitutes a 
reliable test for the presence of an x-ray shutter in 
an ambiguous situation. Strongly decreasing pixel- 

T5 values peripheral from the computed curve there- 
fore indicate that the computed curve does not 
represent an x-ray shutter-edge. 

It is also an object of the invention to provide 
an image processor to process an image formed 

20 by x-ray imaging, which contains a portion cor- 
responding to an x-ray absorption filter device and 
which is arranged to carry out the image process- 
ing method of the invention. This is achieved by an 
image processor according to the invention' which 

25 is characterized in that the processor comprises a 
device for determining an edge of the portion and a 
device for modifying the image by selectively pro- 
cessing the portion. 

Moreover, it is an object of the invention to 

30 provide an x-ray examination apparatus comprising 
an x-ray absorption filter between an x-ray detector 
and an x-ray source for irradiating an object to form 
an x-ray image on the x-ray detector and which is 
arranged to carry out the image processing method 

35 of the invention. This is achieved by an x-ray 
examination apparatus according to the invention, 
characterized in that the apparatus comprises an 
image processor according to Claim 16 and having 
an input port coupled to the x-ray detector. 

40 These and other aspects of the invention will 

become apparent from and will be elucidated with 
reference to the implementations and the embodi- 
ments described hereinafter and with respect to the 
accompanying drawings wherein 

45 Figure 1 shows a diagrammatic representation 
of an implementation of ah image processing 
method in accordance with the invention; 
Figures 2a-c shows an exemplary variation of 
brightness around an x-ray shutter-edge, and 

so first and second spatial derivatives of said vari- 
ation; 

Figure 3 shows a diagram of a division into sub- 
regions in accordance with the image process- 
ing method in accordance with the invention of a 
55 region of interest; 

Figure 4a shows a diagrammatic illustration of 
detection of x-ray shutter-edges in an image 
made by x-ray exposure; 

7 » 
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Rgure 4b shows a result of an image process- 
ing method in accordance with the invention 
having been applied to the image shown in 
Rgure 3b; 

Rgure 5 shows a schematic representation of an 
x-ray examination apparatus comprising a de- 
vice for performing an image processing method 
in accordance with the invention. 
Rgure 6 shows a diagrammatic representation 
of a further implementation of an image pro- 
cessing method in accordance with the inven- 
tion; 

Rgure 1 shows a diagrammatic representation 
of an image processing method in accordance with 
the invention. In block I a region of interest 1 is 
selected^ wherein an x-ray shutter-edge is expect- 
ed. In block II said region of interest 1, is divided 
into sections 2. On the basis of brightness values 
in each of the sections 2, transition-points are de- 
termined, one of them being indicated by the refer- 
ence numeral 3 Furthermore, weights are assigned 
to the transition^points and a curve is fitted through 
the transition-points where use is made of the 
weights for performing the curve fitting so as to 
obtain a fair representation of an x-ray shutter- 
edge. In block IV a fitted curve 60 is shown. Once 
a curve has been determined, further accuracy is 
achieved by computing weights for each of the 
transition-points where for a transition-point a 
weight is computed that depends on the distance 
between the fitted curve and said transition-point. 
Reassigning weights is performed in block V and 
the curve-fitting process is re-iterated under control 
of a comparison step that is represented by block 
VI. When this iteration has converged, reliability of 
the converged curve is computed in block VII on 
the basis of the weight. In this way, because x-ray 
shutter-edges have transition-points that have com- 
paratively high weights whereas anatomical struc- 
tures have transition-points with lower weights, it is 
avoided to a large extent that anatomical structures 
are mistaken for an x-ray shutter-edge. To that end 
a computed reliability is compared in a comparison 
step represented by block VIM with a predeter- 
mined threshold level that is provided from a mem- 
ory means IX. In a final step shown in block X, 
image information of the x-ray shutter-area, being 
the area peripheral to the fitted curve, is deleted. 
The deletion may consist of setting all pixel-values 
on the x-ray shutter-area to a grey-level or a non- 
distracting colour-value that may be either pre- 
determined or can be selected. Thus, said x-ray 
shutter-area is reduced in that image information of 
the x-ray shutter-areas is deleted. As an alternative 
the x-ray shutter-area is reduced in that said x-ray 
shutter-area is discarded from the image. 

Rgures 2a-c shows an exemplary variation of 
brightness around an x-ray shutter-edge, and first 



and second spatial derivatives of said variation. 
Rgure 2a shows dependence of brightness-values 
as a function of position along a line in an image 
made by x-ray exposure. The coordinate x in- 

s creases from the periphery of the image towards 
the centre of the image. As Rgure 2a shows the 
brightness-value B is very low near the periphery of 
the image owing to the presence of an x-ray shut- 
ter. There is a transition to high brightness-values 

10 around the position x t which is representative for 
the edge of an x-ray shutter. Rgure 2b shows the 
first spatial derivative B'(x) in the direction of said 
line of the brightness variations shown in Rgure 2a. 
Notably, at the x-ray shutter-edge x% the derivative 

is B' has a maximum, indicating the brightness vari- 
ations being the most pronounced at the x-ray 
shutter-edge, Rgure 2c shows the second deriva- 
tive B"(x) in the direction of said line. As appears 
from Rgure 2c, the second derivative has a maxi- 

20 mum close to the position x t , where the first deriva- 
tive has a maximum. In relation to the maximum of 
the first derivative, the maximum of the second 
derivative has a position that is towards the periph- 
ery of the image. The feature that the maximum of 

25 the first derivative at the side facing the periphery 
of the image is accompanied by a maximum of the 
second derivative, discriminates a brightness-tran- 
sition having low brightness on the peripheral side 
and having high brightness on the side towards the 

30 centre from a brightness-change having low bright- 
ness on the peripheral side and gradually increas- 
ing brightness towards the centre; such a bright- 
ness-change in practical situations often arises due 
to vignetting in the image. 

35 Rgure 3 shows a diagram of a division into 

sub-regions of a region of interest in accordance 
with the image processing method in according to 
the invention. 

To search for an x-ray shutter-edge of an x-ray 

40 shutter located in the left part of the image first a 
region of interest ROI to tai 1 is selected as is in- 
dicated in Rgure 3. The region ROI to tai extends 
over half the vertical dimension of the image and 
from the left edge of to the middle of the image in 

45 the horizontal direction. In this manner, artifacts in 
the corners are left out of consideration and risk of 
conflict with other x-ray shutters being present is 
diminished. Other x-ray shutters, such as right- 
vertical, top-horizontal or bottom-horizontal can be 

so detected in accordance with the invention by 
choosing regions replacing vertical and horizontal 
and top and bottom. Subsequently, the region 
ROItotai* which consists e.g. of 256x256 pixels is 
divided into elongated strips ROI sub , one of them 

55 being shown of e.g. 256x32 pixels. For each of the 
strips pixel values are averaged in the vertical 
direction, i.e. to each of the 256 positions in the 
longitudinal direction of the strip currently being 
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processed the average of the 32 pixel-values hav- 
ing the same horizontaJ position is assigned. This 
yields a dependence of (average) brightness along 
the elongated direction of the current strip, as is 
depicted e.g. in Figure 2a. 

Figure 4a shows a diagrammatic illustration of 
detection of x-ray shutter-edges in an image made 
by x-ray exposure. Figure 4a shows ah image 
made by x-ray exposure featuring various anatomi- 
cal structures such as a bone 10 and a artery 11 
and a catheter 12. Furthermore, several elongated 
strips 21, 22 and 23, similar to the strips 2 in 
Figure 3 are indicated. In each of the strips several 
brightness transitions occur. These brightness-tran- 
sitions originate from either anatomical structures 
: or from an x-ray shutter. E.g. in strip 21, at position 
31 brightness changes from high to low when pro- 
ceeding towards the periphery of the image, at 
position 32 the brightness B increases to a higher 
value and finally at position 33 brightness again 
drops to a very low value. Consequently, both 
brightness transitions at 31 and 33 are accom- 
panied by respective maxima of the spatial deriva- 
tives B' and 8", as discussed with reference to 
Figures 2a-c. On the other hand the brightness 
transition at position 32 is accompanied by a mini- 
mum of B'; so that only the transitions at 31 and 
33, respectively are considered candidates for a 
sought x-ray shutter-position. Because the transi- 
tion at position 31 is accompanied by another 
transition being more peripheral, viz. the transition 
at position 33. to position 31 a smaller weight-factor 
is assigned and to the position 33 a large weight- 
factor is assigned. 

Similarly, two further strips 22 and 23 are in- 
dicated, each of them comprising several bright- 
ness transitions at positions 41, 42, 43 and 51, 52, 
53, respectively. Of these, the transitions at 43 and 
53 are assigned a large weight-factor value. By a 
known curve-fitting method a curve 60 is drawn 
through the group of transition-points formed by 
(31, 33, 41, 43, 51. 53). The curve 60 forms a fair 
estimate of the x-ray shutter-edge of a left x-ray 
shutter. 

The accuracy of the fitted curve which is repre- 
sentative for an x-ray shutter-edge is increased by 
iterating. After an estimate is obtained, weight-fac- 
tors are reassigned to the positions that were used 
for obtaining the current fitted curve. The process 
is stared by determining peak-points in each of the 
sections and choosing a curve that connects said 
peak-points as a first estimate for a curve repre- 
senting an x-ray shutter-edge at issue. A peak- 
point is a transition-point such that the average 
distance between a predetermined number of 
neighbouring transition-points and said peak-point 
is the smallest; empirically it has appeared that the 
choice of five neighbouring transition-points is very 



adequate to determine a peak-point. When' re-as- 
signing weight-factors, distances between said por 
sitions and the current fitted curve are taken into 
account such that positions that are close to the 

s current fitted curve are given higher weight-factors 
than position that are more distant from the current 
fitted curve, for example by means of a Gaussian 
weight function. Then, on the basis of new values 
of weight-factors being assigned, again positions 

10 being most probable candidates for lying on an x- 
ray shutter-edge are selected and curve fitting is 
carried-out employing said new weight-factor val- 
ues. 

Still further refinement of the reliability of the 
75 image processing method in accordance with the 
invention is achieved by considering image fea- 
tures around a transition-point, i.e. a presumed x- 
ray shutter-edge. A first feature parameter of a 
brightness-transition is the height of the transition 
20 which is determined by the difference between an 
averaged brightness values at the side of the tran- 
sition towards the periphery of the image and an 
averaged brightness at the side of the transition 
towards the centre of the image. The averaging is 
25 e.g. performed over a predetermined number of 
pixels at a relevant side of a transition at issue. 

A second feature parameter is the roughness 
of the image at the peripheral side of a brightness- 
transition at issue. In practical situations, x-ray 
30 shutters appear distinct from anatomical structures 
in that within an x^ray shutter-area there is very 
little structure; this is quantified by image rough- 
ness which is defined as the statistical variance of 
brightness-values in an area of the image af issue. 
35 An attractive alternative for quantifying roughness 
is obtained on the basis of second spatial deriva- 
tives of brightness which are interpreted as the 
difference between a local brightness-value and a 
running average built-up by three neighbouring 
40 pixel-positions in total. Then roughness is defined 
quantitatively as a reciprocal of an average over 
the modulus of second spatial derivatives in an 
area of the image at issue. Low roughness values 
then correspond to an area of the image having 
45 little image structure, thus having a high probability 
of being an x-ray shutter-area. Further an image 
processing method in accordance with the inven- 
tion is now refined in that weight-factors that are 
assigned to positions that are recognised as can- 
so didates for being on an x-ray shutter-edge are 
made increasing as a function of the height of the 
transition and are made decreasing as a function of 
the roughness in the area of the image peripheral 
of the candidate at issue. 
55 After all x-ray shutter-edges in the image have 

been determined, image-information in the areas of 
the image peripheral to any of the curves cor- 
responding to x-ray shutter-edges is deleted from 
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the image and as a result a processed image 
comprising only relevant medical information Is ob- 
tained. In particular, in areas corresponding to x-ray 
shutters as determined according to the invention, 
all pixels are set to the same value; in this manner 
x-ray shutter-areas are blanked. Figure 4b shows a 
processed image obtained from the image shown 
in Figure 4b. 

The processed image has as an advantage that 
it does not include parts that are distracting when 
examining the image. Furthermore, when a hard- 
copy e.g. on film of a plurality of images having 
been processed according to the invention is 
made, an advantage is achieved in that more eco- 
nomic use is made of hard-copy material. 

Figure 5 shows a schematic representation of 
an x-ray examination apparatus comprising a de- 
vice for performing the image processing method 
in accordance with the invention. An x-ray, source 1 
is provided for irradiating an object 2 so that an 
image carrying x-ray beam is formed. By means of 
an x-ray detector 4 and an image read-out device 5 
an image made by x-ray exposure is formed an 
supplied to an image-memory 6. An image proces- 
sor 7 is employed for converting the image stored 
in the image-memory 6 into a processed image 
from which x-ray shutter-areas are deleted. The 
processed image is supplied to an image buffer- 
circuit 8 for further processing such as display on a 
monitor screen, formation of a hard-copy or to be 
archived. 

The image processor 7 will now be described 
in more detail. A selector 11 selects a region of 
interest ROI tot ai wherein an x^ray shutter-edge is to 
be detected. By a divider 12 the region of interest 
is divided into elongated strips ROIsu b . The selector 
11 and the divider 12 are controlled by an iterator 
14 for subsequently selecting said elongated strips 
for determining an x-ray shutter-edge and for sub- 
sequently selecting regions of interest for determin- 
ing various x-ray shutter-edges. Brightness-values 
of pixels in an elongated strip are averaged along 
the transverse direction of the elongated strip at 
issue so that a one-dimensional brightness distribu- 
tion being representative for image information 
along the longitudinal direction of the elongated 
strip, is formed by an averager 13. First and , sec- 
ond derivatives of the pne-dimensional brightness 
distribution are computed by differentiators 15 and 
16, respectively. From the output of the differen- 
tiators candidate-positions for x-ray shutter-edge 
positions are determined by a comparator 17. To 
candidate x-ray shutter-edge positions weight-fac- 
tors are assigned by way of a weighting-means 18. 
Said weight-factors are computed by the weighting 
means 18. To that end the height of a brightness- 
transition at a positions at issue is computed by 
subtracter 19; the roughness associated to the 



brightness-transition at issue is computed by an 
evaluator 20. Through the candidate-x-ray shutter 
positions in successive elongated strips and using 
their associated weight-factors candidate x-ray 

5 shutter positions are selected and a curve is fitted 
by a curve-fitter 21 . In addition to producing a 
curve that fits the selected positions, the curve- 
fitter computes updated weight-factors depending 
on the distances between respective selected posi- 

io tions and the fitted curve. Under control of the 
iterator 14, selection of candidate x-ray shutter- 
positions and subsequent curve-fitting is repeated 
so as to obtain a more reliable result for a curve 
being a fair approximation of an x-ray shutter^edge. 

/5 The curve found after iteration has converged is 
supplied to a verification-means 22 for to verify that 
said curve is indeed representative of an x-ray 
shutter-edge. The verification-means comprises a 
computation unit for calculating an average dis- 

20 tance to a curve at issue of a predetermined num- 
ber of transition-points. In particular five transition- 
points are chosen that have the smallest distances 
to the curve at issue. Moreover, the verification 
means comprises detection means for ascertaining 

25 that for a curve at issue, an x-ray shutter-edge is 
present in an opposite region of interest in the 
image. Subsequent to verification of approximating 
curves for various x-ray shutter-edges being 
present in the image stored in the image-memory 6 

30 , parts of the image being peripheral to any said x- 
ray shutter edges are deleted by an eraser 23 and 
the processed image obtained in this manner is 
supplied to the image buffer-circuit 8. 

It is noted that in an x-ray examination appara- 

35 tus in accordance with the invention, the functions 
of the image processor may be performed by a 
suitably programmed computer or by a special 
purpose processor having circuit-means that are 
arranged to perform said functions. 

40 Figure 6 shows a diagrammatic representation 

of a further implementation of an image processing 
method in accordance with the invention. In block 
A a rough detection of transition-points is per- 
formed. In this phase of the image processing 

45 method the full image having, for example, 
512x512 pixels is subsampled down, for example, 
to 64x64 pixels by taking one point out of eight in 
the horizontal arid vertical directions. This consider- 
ably reduces the computational burden at the cost 

so of a lesser accuracy in localisation of the transition- 
points. Accuracy will be recovered in the sequel of 
i the image processing method, viz. in block B. The 
following computational steps are performed. The 
subsampled image is smoothed using a straight 

55 vertical kernel. This is done by adding to each 
pixel-value the pixel-value of pixels in the neigh- 
bourhood in said sub-sampled imaged In this way, 
sharp horizontal transitions are retained in this way 
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while vertical transitions are smoothed out. 

Subsequently, nine rows are selected as sec- 
tions, out of the 64 rows of the subsampled image: 
A systematic search for transition-points that lie on 
an x-ray shutter-edge is then performed. To that 
end, in each row the average pixel-value of pixels 
located between a transition-point and a periphery 
of the image is estimated by taking the sum of 
intensities of four pixels located between the pixel 
for which said average is to be computed and the 
periphery of the image. Next, the component in a 
direction orthogonal to a periphery of the image of 
a gradient of pixel-values is evaluated by comput- 
ing differences between consecutive points in said 
section. 

In block B accurate localisation of transition- 
points is represented. From the phase performed in 
block A, nine roughly localised transition-points are 
provided. The phase of block B aims at increasing 
accuracy of localisation of transition-points in the 
512x512 image. A local search around each of 
said nine transition-points is performed within the 
same section for a neighbourhood consisting of, for 
example, 17 pixels, viz. eight on either side of a 
transition-point at issue. The local search com- 
prises smoothing in a direction parallel to said 
periphery of the image by adding pixel-values that 
neighbour the transition-point at issue in a direction 
parallel to said periphery* Subsequently, the sum is 
computed of four pixel-values of four pixels located 
between said periphery and the transition-point at 
issue, said sum being representative for an aver- 
age brightness towards said periphery. A gradient 
in a direction orthogonal to said periphery is com- 
puted by calculating differences of pixel-values of 
the transition-point at issue and two of its neigh- 
bours that are located between said transition-point 
and the centre of the image. Next, the ratio of said 
gradient and said average is computed for each 
transition-point in each of the sections. Accurate 
transition-points having a comparatively high prob- 
ability for being on as x-ray shutter-edge are se- 
lected as those transition-points having a value for 
said ratio being the highest in said neighbourhood 
of 17 pixels. 

Block C represents the fitting of a curve, nota- 
bly the best fit straight line through the accurate 
transition-points that are provided from the phase 
of block B. Accurate transition-points which are 
comparatively far from said straight line are consid- 
ered as outliers and are rejected. Alteratively, tran- 
sition points can be weighted according to the 
distance to the curve, for example by a Gaussian 
weighting function. The computational procedure 
for the phase of block C is iterative; it represents a 
compromise between accuracy and computation 
time optimisation. The estimation of the X-ray shut- 
ter edge position in the image and angle with a 



predetermined axis of the image is then refined by 
determining the best fit straight line to the most 
likely accurate transition-points. For the following of 
the discussion, said predetermined axis will be 

6 chosen to be a vertical axis of the image, although 
the implementation is not at all restricted to this 
choice. The pixels of the image are arranged as a 
square matrix of columns and rows. The x-ray 
shutter-edge is first assumed to be vertical, its 

10 column number is taken to be that of the accurate 
- transition-point with the largest number of substan- 
tially vertically aligned accurate transition-points. 
Each of the accurate transition-points is assigned a 
likelihood score. Let Ax represent the maximum 

is change in column number of a said accurate transi- 
tion-points to account for a maximum angle depar- 
ture of, for example, ±5 • from a vertical edge. The 
likelihood score of any one of the accurate transi- 
tion-points is calculated by adding up the number 

20 of accurate transition-points with column numbers 
differing by less than Ax plus those with column 
number differing by less than Ax/2 plus those with 
column numbers differing by less than Ax/4. The 
first guess for the position of the x-ray shutter-edge 

25 is taken to be a vertical line at the column number 
of the accurate transition-point having the highest 
likelihood score. In practice Ax is chosen to be a 
width of eight pixels. 

To proceed further -with a low computational 

30 burden and to decrease required computation time, 
accurate transition-points which lie comparatively 
far away from said guessed x-ray shutter edge 
position are rejected, for example according to the 
following test-rules. If four or more accurate transi- 

35 tion-points are within a distance of Ax/4 from said 
guessed x-ray shutter-edge position then accurate 
transition-points having a distance larger than Ax/2 
are rejected. If this test-rule does not . apply, the 
test-rule is repeated after replacing the value of Ax 

40 by 2Ax. If neither of these test-rules apply, i.e. less 
than four accurate transition-points within a dis- 
tance of Ax from said guessed x-ray shutter-edge 
position, it is concluded that no x-ray shutter is 
present. 

45 The accurate transition-points that have been 

retained after applying said test-rules are employed 
to more accurately localise the x-ray shutter-edge 
position in an iterative way under the control of an 
iteration control indicated as block E. The initial 

50 guess of the x-ray shutter-edge is a vertical line at 
the column number having the highest likelihood. 
/For each of the retained accurate transitionrpoints 
respective normal distances to said vertical line are 
computed and Gaussian weights depending on 

55 said normal distances are assigned such that the 
extreme points for an x-ray shutter having a tilt of 
5* with the vertical axis are given weights of i. 
The function of these weights is to diminish the 
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influence of those retained accurate transition- 
points which are too poorly aligned with respect to 
others. A conventional least mean square regres- 
sion is employed to determine a line that is repre- 
sentative for a new estimate of an x-ray shutter- 
edge position. In the iteration process, normal dis- 
tances of retained accurate transition-points are 
recalculated for the new estimate for the x-ray 
shutter-edge position. Subsequently, new Gaussian 
weights are assigned where a Gaussian distribution 
with a smaller width is used. Following that, a new 
best fit straight line is determined, thus improving 
the estimate for the position of the x-ray shutter 
edge. 

After a line that represents an x-ray shutter 
edge position has been provided by the phase of 
blocks C and E by iteration a decision-taking phase 
follows as is represented by block D. After the 
iteration has converged the root mean square dis- 
tance Ap of the retained accurate transition-points 
is calculated. The following situation may then oc- 
cur. If Ap is larger than 2.6 pixels it is decided that 
no x-ray shutter is present. The image processing 
method is then terminated. The value of 2.6 pixels 
has been determined empirically and may be dif- 
ferent in other implementations or applications. If 
all nine accurate transition-points are retained and 
Ap is less than 1.2 pixels, it is decided that there is 
an x-ray shutter and the image processing method 
proceeds by reducing the area 61 that is peripheral 
to the line provided by the phase of block D. If only 
eight accurate transition-points are retained and Ap 
is less than 0.8 pixels it is also decided that there 
is an x-ray shutter and subsequently the x-ray 
shutter-area is reduced in block E. If none of the 
above mentioned situations occur the presence of 
an x-ray shutter is ambiguous and three further 
tests are performed In block D. 

If such an ambiguous situations occurs, there 
are provided from previous phases a set of ac- 
curate transition-points. To save computation time 
further tests are carried-out using the 64x64 sub- 
sampled image. The mean pixel-value along two 
vertical line segments an inner one lying towards 
the centre of the image and the other, outer one, 
lying towards said periphery are computed, the 
common length of these two line segments is taken 
to be slightly smaller than an estimated length of 
an x-ray shutter and to be compatible with the 
rounded shape of the x-ray image intensifier. A 
high mean pixel-value on the inner-segment and a 
low value on the outer-segment is compatible with 
the presence of an x-ray shutter. The quantity 
computed in this test is the ratio R rt of said two 
mean pixel-values and that quantity is compared 
with a predetermined reference value to decide on , 
the presence of an x-ray shutter. 



As second test, the maximum value L max of the 
second spatial derivative (Laplacian) of the pixel- 
values on the outer segment If this maximum 
value is found to be large, it indicates the presence 

5 of a anatomical structure towards the periphery of 
the accurate transition-point at issue and hence the 
likelihood of the presence of an x-ray shutter is 
reduced. A reliable criterion is a linear combination 
M of the ratio R rt and the maximum value L max . The 

io coefficients involved in this liner combination are 
those optimising the statistical recognition results 
on the part of the image used for training the 
algorithm. The linear combination M is compared to 
a predetermined reference value so as to decide 

75 on the presence of an x-ray shutter. 

A final independent test that is used in an 
ambiguous situation also aims at rejecting false x- 
ray shutter detections. In this final independent test 
the maximum negative slope on a line orthogonal 

20 to said periphery of an image is computed between 
an accurate transition-point at issue and said pe- 
riphery. If this maximum negative slope is larger 
than an predetermined reference value it is de- 
cided that no x-ray shutter is present. As pixel 

25 values hardy ever decrease strongly in a direction 
of an x-ray shutter, considering said maximum neg- 
ative slope constitutes a reliable test for the pres- 
ence of an x-ray shutter in an ambiguous situation. 
Strongly decreasing pixel-values peripheral from 

30 the computed curve therefore indicate that the 
computed curve does not represent an x-ray shut- 
ter-edge. 

Claims 

35 

1. An image processing method to process an 
image formed by x-ray imaging, the image 
containing a portion corresponding to an x-ray 
absorption filter, characterized in that the im- 

40 age processing method comprises 

- modifying the image by selectively pro- 
cessing said portion. 

2. A image processing method as claimed in 
45 Claim 1, characterized in that the processing 

comprises 

- determining an edge of the portion for 
selecting said portion as a part of the 
image between the edge and a periphery 

so of the imaige. 

3. An image processing method as claimed in 
Claim 1 or 2, characterized in that the process- 
ing comprises: 

55 - removing at least a part of said portion 

from the image. 
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4. An image processing method as claimed in 
Claim 1, characterized in that the processing 
comprises: 

- assigning substantially uniform values to 
pixels in the portion in the image. s 

5. An image processing method as claimed in 
any one of Claim 1 or Claim 2, characterised in 
that the detection of the edge comprises 

a grouping step to select a group of transition- 10 
points in the image having a higher brightness 
towards the centre of the image and a lower 
brightness towards the periphery of the image, 
and 

a curve-fitting step to fit a curve through said 75 
group of transitioo-points, said curve forming a 
representation of said edge. 

6. An image processing method as claimed in 
Claim 5, characterised in that the grouping 20 
step comprises 

a region-selection step to select a region of 
interest extending from a periphery of the im- 
age and 

a division-step to divide said region of interest 25 
into sections and 

a weight-assignment step to assign weights to 
said transition-points, a weight for a transition- 
point depending on image features around said 
transition-point and a weight for a transition- so 
point depending decreasingly on a distance 
between said transition-point and said periph- 
ery. 

7. An image processing method as claimed in 35 
Claim 5 or 6, characterised in that the grouping 
step comprises a first-estimate step to deter- 
mine a first-estimate for a position of a transi- 
tion-point in a section and an accurating step 

to determine a transitiorv^point from said first- 40 
estimate. 

8. An image processing method as claimed in 
Claim 6 or 7, characterised in that the weight- 
assignment step comprises 45 
computation of image roughness, and 
computation of transition-height, and 
computation of said weights depending on im- 
age roughness and transition height for respec- 
tive transition-points. so 

9. An image processing method a£ claimed in 
any one of Claims 5 to 8, characterised in that 
the weight-assignment step comprises 
computation of first and second spatial derjva- 55 
tives of brightness-values near respective tran- 
sition-points and 

employing substantial concurrence of a maxi- 



mum of said first derivative with a brightness- 
variation and with a maximum of said second 
derivative at the peripheral side of said bright- 
ness^variation for selecting the position of said 
brightness-variation as a transition-point and 
assigning said weight on the basis of said 
concurrence. 

10. An image processing method as claimed in 
any one of Claims 5 to 8, characterised in that 
the weight-assignment step comprises 
computation of first spatial derivatives of 
brightness-values peripheral to respective tran- 
sition-points, computation of respective aver- 
age brightness-value of pixels towards the cen- 
tre to said transition-points, computation of re- 
spective ratios of said respective first spatial 
derivatives and said respective average bright- 
ness-values and assigning said weights on the 
basis of said ratios. 

11. An image processing method as claimed in 
any one of Claims 5 to 10, characterised in 
that 

the weight assignment step is followed by a 
point-selection step for selecting a transition- 
point having the highest weight in a section 
and retaining only said selected transition- 
points in said selected group. 

12. An image processing method as claimed in 
any one of Claims 5 to 11. characterised in 
that 

the curve-fitting step is followed by a reassig- 
ning step for reassigning a value to a weight of 
a transition-point in dependence on a distance 
between said transition-point and a curve fitted 
through the group whereto said transition-point 
belongs, and in that the curve-fitting step is re- 
iterated employing said re-assigned weight. 

13. An image processing method as claimed in 
any one of Claims 2 to 12, characterised in 
that determining the edge is followed by an 
edge-verification step comprising computation 
of a reliability-parameter for an edge being 
determined and comparison of said reliability- 
parameter with a predetermined threshold lev- 
el. 

14. An image processing method as claimed in 
Claim 13, characterized in that for a plurality of 
transition-points said reliability parameter is 
computed and in a situation where comparison 
of said reliability parameter to said, predeter- 
mined reference level indicates the presence 
of said portion corresponding to an x-ray ab- 
sorption filter for less transition-points than a 
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predetermined reference number the edge-ver- 
ification step is supplemented by comparison 
of pixel-values of collections of pixels on dif- 
ferent sides of the curve. 

15. An image processing method as claimed in 
Claim 13, characterized in that plurality of tran- 
sition-points said reliability parameter is com- 
puted and in a situation where comparison of 
said reliability parameter to said predetermined 10 
reference level indicates the presence of said 
portion corresponding to an x-ray absorption 
filter for less transition-points than a predeter- 
mined reference number the edge-verification 
step is supplemented by determining a mag- is 
nitude of a first spatial derivatives of bright- 
ness-values of pixels in the image between the 
curve and the periphery of the image, said 
derivative being negative in the direction from 

the curve to the periphery of the image and 20 
comparing said magnitude to a predetermined 
reference value. 

16. An image processor to process an image 
formed by x-ray imaging, which contains a 25 
portion corresponding to an x-ray absorption 
filter, characterized in that the processor com- 
prises a device for 

. - determining an edge of the portion 
- and a device for modifying the image by 30 
selectively processing the portion. 

17. An x-ray examination apparatus comprising an 
x-ray absorption filter between an x-ray detec- 
tor and an x-ray source for irradiating an object 35 
to form an x-ray image on the x-ray detector, 
characterized in that the apparatus comprises 

an image processor according to Claim 16 and 
having an input port coupled to the x-ray de- 
tector. 40 
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